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ABSTRACT 

We study the alignments between the angular momentum of individual objects and 
Q ■ the large-scale structure in cosmological numerical simulations and real data from the 

Sloan Digital Sky Survey, Data Release 6. To this end we measure anisotropies in the 
C/3 , two point cross-correlation function around simulated halos and observed galaxies, 

studying separately the 1- and 2-halo regimes. The alignment of the angular momen- 
tum of dark-matter haloes in ACDM simulations is found to be dependent on scale 
^SJ ■ and halo mass. At large distances (2-halo regime), the spins of high mass haloes are 

^ ' preferentially oriented in the direction perpendicular to the distribution of matter; 

\ lower mass systems show a weaker trend that may even reverse to show an angular 

■ momentum in the plane of the matter distribution. In the 1-halo term regime, the an- 
' gular momentum is aligned in the direction perpendicular to the matter distribution; 
, the effect is stronger than for the 1-halo term and increases for higher mass systems. 

' On the observational side, we focus our study on galaxies in the Sloan Digital Sky 

' Survey, Data Release 6 (SDSS-DR6) with elongated apparent shapes, and study align- 

OO , ments with respect to the major semi-axis. We study five samples of edge-on galaxies, 

■ the full SDSS-DR6 edge-on sample, bright galaxies, faint galaxies, red galaxies, and 

blue galaxies (the latter two consisting mainly of ellipticals and spirals respectively). 
Using the 2-halo term of the projected correlation function, we find an excess of struc- 

' ture in the direction of the major semi-axis for all samples; the red sample shows the 

H \ highest alignment (2.7 ±0.8%) and indicates that the angular momentum of flattened 

spheroidals tends to be perpendicular to the large-scale structure. These results are 
in qualitative agreement with the numerical simulation results indicating that the an- 
gular momentum of galaxies could be built up as in the Tidal Torque scenario. The 
1-halo term only shows a significant alignment for blue spirals (1.0 ± 0.4%), consis- 
tent with the 1-halo results from the simulation but with a lower amplitude. This 
could indicate that even though the structure traced by galaxies is adequate to study 
large-scale structure alignments, this would not be the case for the inner structure of 
low mass haloes, M ^ lO^^/i^^M©, an effect apparently more important around red 
g — r > 0.7 galaxies. 
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1 INTRODUCTION 

In the current paradigm of large-scale structure formation 
in the Universe it is well established that the angular mo- 
mentum of dark-matter (DM) halos is determined early in 
the history of the Universe by the gravitational interac- 
tion between the quadrupole of the collapsing proto-halo 



region and the surrounding matter. This was first formu- 
lated fo r the h ierarchical the o ry of structure formation by 
,Whit"3 (|1984| ). iDoroshkeviciil (|l970l ). and iPeeblesi (j 19691 ) . 
and it is claimed that this is a natural consequence of the 
pertur bative treatment of the gravitation al instability sce- 
nario (jPorciani. Dekel fc HoflmanI l|2002cj )). Several inher- 
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ent properties to the tidal torque process (Tidal Torque 
Theory, TTT) have been noticed and studied in the lit- 
erature. For instance it is expected from TTT, and also 
supported to some degree by N-body simulations, that 
the angular momentum s t arts growinR linearly with time 
ijCatelan fc Theunslll996bl : |Peeblesl[l969l : IWhitelll984h . and 
that later on this process looses efficiency after the halo turn- 
around. This is due to the reduction of halo inertia by col- 
lapse, and to the continued expansion of the neighboring 
matter responsible for the tidal torque. 

Th ere are sever al analytic studies o f the TTT sce- 
nario llWhitd 11984 iDoroshkevichI Il970l : iPeeblesI 1 19691 : 
ICatelan fc Porcianil l200ll), a s well as semi-analytic ap- 
proaches ( Catelan fc TheunsI Il996ah . and fully non-linear 
work using cosmolo gical simulations (iBarnes fc Efstathiou 



our Galaxy, and more r ecently on large galaxy surveys 
llHirata et atl 12001 l2004l: iLee fc Erdogdull2007l: iLee fc PenI 



1987 



200C 



Porciani j^ Dekel fc Hoffmanj 



Peirani. Mohavaee. fc de Freitas Pachecol 



12004 



Sugerman. Summer s fc Kamionkowski 

2002allbl : 



Nagashima fc Goudal 19981 ). The agreement is that the 



resulting halo population is poorly rotationally supported; 
the velocity dispersion at virialisation is found to be the 



most important mechanism for halo equilibrium. As a result, 
the distribution of spin parameter^J i s log-normal which 
peaks at quite l o w values, Amed — 0.03 5 (iBullock et al.ll200ll: 



,t quite l ( , 

Gardne^ |200ll: iD'Onghia fc Burkertl l2004l : iMaccio' et all 



20061 : iBett et al.l l2007<i. Unfortunately it is not possible to 



obtain a direct meas urement of spin parame ters in actual 
galaxies. Despite this, [Hernandez et al] (|2007l ) used indirect 
methods to derive distributions of galaxy spin parameters 
consistent with theoretical predictions. 

From an observational point of view, it can be easier 
in some cases to determine the direction of the angular 
momentum of a given galaxy rather than its amplitude. 
There has been a recent increase in the interest concern- 
ing the possibility of galaxy alignments with the large 
scale structure. Such studies are important to test the 
already well developed TTT, and to asses the possibility 
that such correlations affect the results of cosmic shear 
measurements. Regarding the former, TTT produces noisy 
predictions due to the non-linear nature of the accre- 
tion processes of clumpy matter (|Vitvitska et al.l |2002| : 
iPorciani. Dekel fc HoffmanI l2002ah . on scales below the 
coherent scale of flows predicted by linear theory. The 
second tests for the effects on cosmic shear studies used 
to analyse weak lensing, may answer whether there is an 
intrinsic alignment in the structure traced by galaxies 
that could result i n a spurious c ontrib u tion to the shear 
power spectrum llHirata et all l2007l: iBridle fc Abdalla 



2007'; Takada fc White"2004'; 'King , 
2005; Hcymans fc Heavens 2003; 



Schneider^OO 
Hcymans et al 



lT 



.King 
2004 



Brown et all |2002|: [ Catelan. Kamionkowski fc BlandfordI 



20011 : ICrittenden et all |2001| : ICroft fc Metzled l2000l : 
Jind |2002| ). Several authors have searc hed for galaxy 



align ments in the super-galactic plane (|de Vaucouleurd 
I1953I ). extending out to several tens of megaparsecs from 



^ The dimensionless spin parameter as defined bv lPeebleS (Il969l ) 
is A = L^\E\/GM^/'^ , where E is the total internal energy of 
the halo and L and M are its angular momentum and mass re- 
spectively. This number can be understood as a measure of the 
importance of rotational velocity vs. velocity dispersion as the 
dominant support mechanism of a halo 



I2OO7I : iMandelbaum et al.l I2OO6I ). iLee fc Erdogdul (|2007D 
claim a direct evidence of an alignment between galaxy 
position angle in 2MRs and the tidal field deduced from 
galaxy positions. 

An interesting approach to the analysis of galaxy align- 
ments is to use the radial direction to the centres of voids 
as characteristic directions in the large scale structure. It 
is well known that as they evolve, large voids tend to be 
rounder, a behaviour which is opposite to that shown by 
bound structure. Negative density perturbations in the ini- 
tial fluctuation field that will later become voids, are char- 
acterised by a decreasing aspheri city as the expansion of the 
Universe proceeds, as is shown bv lSheth fc van de WevgaertI 
l|2004l ) (and references therein). Taking advantage of this 
particular characteristic of voids, it is possible to study the 
whether galaxy or halo angular momenta are aligned with 
the void shells just by considering the angle with respect to 
th e direction to the void centr e . This is the approach used 
bv lTruiillo. Carretero fc Patiril (|2006l ) applie d to the SDSS 
and 2 d FGRS galaxy r edshif t sur veys, and bylBrunino et al.l 
(|2007l ). [Cuesta et al l (l2008li andlPatiri et a\U2QQ^ t o nu- 
merical simulations. iTruiillo. Carretero fc Patiril ( 20061 ) find 
with a high level of confidence that spiral galaxies located 
on the shells of the largest cosmic voids have rotation 
axes that lie preferentially on the void surface, in agree - 
ment with the simu l ation results of Brunino et al.l |20o3); 
ICuesta et al] (j2008l ): IPatiri et af] (|2006l ). Other more gen- 
eral methods that also take into account alignments with 
large sc ale filaments in numerical simulations show similar 
resuhs (jAragon-Calvo et aLllioOTl : iHahn et al.ll2007al lbh. 

In this paper we present results from a study of align- 
ments between individual objects and their surrounding 
structure in quantitative concordance with TTT predictions 
and previous observational results. We use a novel approach 
to study alignments with the large-scale structure using as 
a reference the galaxy /dark matter halo angular momen- 
tum. Using numerical simulations we study alignments pre- 
dictions in a Cold Dark Matter model with a cosmological 
constant (ACDM), and compare such predictions with the 
results obtained from the Sloan Digital Sky Survey, Data 
Release 6 (SDSS-DR6, Adelman-McCarthy et al., 2007). To 
this end we use the galaxy two point correlation function 
applied both, in its three-dimensional definition to the full 
numerical simulation and in its projected version to mimic 
the observational limitations. The latter allows us to de- 
termine that we should be able to detect the expected low 
signal alignments, even in small galaxy samples. This allows 
us to measure the alignment of galaxies and halos over wide 
ranges in mass and separation, and also to quantify the de- 
pendence of this alignment on halo mass, galaxy colour and 
luminosity. 

This paper is organised as follows. Section 2 presents 
the numerical simulation and the statistical tools we will 
use in this work; it also presents an analysis of the impact 
of the observational biases on the detection of alignments. 
Section 3 contains the analysis of observational samples of 
galaxies from the Sloan Digital Sky Survey. The comparison 
between models and observations as well as the discussion of 
the results is presented in Section 4, and Section 5 contains 
the final conclusions from this work. 
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2 ANISOTROPIES IN THE 3D CORRELATION 
FUNCTION 

In this section we analyse the three-dimensional two-point 
correlation function, using cluster-particle pairs, measured 
using ACDM numerical simulations, taking into account the 
orientation with respect to the centre halo angular momen- 
tum. To this end, we have run two periodic simulations of 
500 and 60 /i~^Mpc computational box sizes, for a flat low- 
density Universe, with a matter density f2m = 1 — = 0.3, 
Hubble constant Ho — 70 km s~^ Mpc~^, and normal- 
isation parameter as — 0.8. The particle resolutions are 
rup J5 7.2 X 10^" Mq for the large simulation, and 
m-p 1.2 X 10* Mq for the small one. The identification 
of particle clumps has been carried out by means of a stan- 
dard friends-of-friends algorithm with a percolation length 
given by / = 0.17 u'^'^, where V is the mean number den- 
sity of DM particles. Both simulations have been performed 
using the first version of the GADGET code developed by 
Springel et al. (2001). We consider halos with at least 30 
and 100 particles in the large and small simulations, respec- 
tively, in order to ensure an accurate measurement of the 
halo angular momentum. 

The spatial halo-particle cross-correlation function ^(r) 
which we will use in this section, measures the excess prob- 
ability with respect to a random distribution, that a DM 
particle will reside at a distance r away from a given halo 
centre, in the volume element dV , 
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Figure 2. Ratios between correlation functions for neighbors 
in the directions perpendicular and parallel to the halo angular 
momentum. Diamonds (triangles) correspond to the ratio over 
the 2-halo (1-halo) term scales. 



dP = i}[l + i{r)] dV. 



(1) 



A standard method to measure ^(r) in a simulation box is 
to count pairs at a given distance bin, and then normalise 
this number by the expected number of pairs for a constant 
number density of objects. In our measurements we perform 
these counts stacking pairs depending on the angle between 
their position relative to the centre halo and its angular mo- 
mentum. We consider tree cases, (i) using all pairs regardless 
of this angle (spherical shell volumes, isotropic case), (ii) us- 
ing all pairs subtending an angle from the angular momen- 
tum less than a given threshold 9i (parallel case), and (iii) 
adding pairs at lower inclinations than a limit angle 62 from 
the perpendicular plane to the angular momentum (perpen- 
dicular case) . We choose the threshold angles, for the parallel 
and perpendicular cases mentioned above {9i_ and 02, respec- 
tively), so that the volumes in each case are the same. This 
can be achieved by setting, sin(S2) = 1 — cos(0i) = x, and 
choosing a value for the threshold parameter x- Selecting 
X ^ 0.5 implies angles Ox sC 60° and 62 ^ 30°. Throughout 
this paper we use x ~ 0-5. 

Errors in the correlation function are estimated using 
the jacknife technique, with a total of 50 subsamples for the 
numerical simulation (20 jacknifes for the SDSS data anal- 
ysis). We find that our errors are stable for this number of 
jacknife subsamples (we tested from 20 to 100) and does 
not produce underestimations. In the case of quotients be- 
tween correlation functions, we apply the jacknife technique 
to the actual measurement of the ratio. We adopt the jack- 
nife technique since it has been reported to provide similar 
error estimates to those obtained from independent samples 
(see for instance, Padilla et al., 2004). 



2.1 Simulation results 

The analysis of alignments of structure with the angular 
momentum of DM haloes presents a first difficulty in that 
in order to obtain a reliable and stable direction for the 
halo spin, it is often necessary to use a large number of 
particles per halo. We find that using 30 particles provides 
reasonable results, but this imposes a rather high lower limit 
to halo masses of 2.2 x lO^'^ Mq for the large simulation, 
and 3.6 x 10^ h'^ Mq for the smaller one. We analised the 
distribution of angular momentum parameters and find that 
the reported high A tail (Bett et al., 2007) due to low number 
of particles per halo, is not present in any of our samples of 
DM haloes. We notice that the clustering signal in the small 
simulation on large scales may be afi'ected by the small box 
size. 

From this point on, we will centre most of our analysis 
on samples of DM haloes with different masses. From the 
small simulation we draw four DM halo subsamples sepa- 
rating haloes of different mass starting at a minimum of 100 
particles or M ^ 1.2 x 10^" Mq (the lowest mass sam- 
ple considered includes haloes of up to 300 particles each; we 
do not use haloes of lower masses even though we are still 
above the minimum number of particles we consider safe to 
consider). Nine samples are drawn from the large simula- 
tion and contain DM haloes with more than 30 particles. 
The subsample with the most massive haloes starts at 1000 
particle members corresponding to M ^ 7.2 x 10^'' Mq. 

The left panels of Figure [1] show the cross-correlation 
function between low mass DM haloes and DM particles 
(top-left) in the directions parallel and perpendicular to 
their angular momenta, and in all directions (diamonds. 
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Figure 1. Left panel: Spatial Halo— Dark-matter correlation function of low mass haloes in the small numerical simulation. The triangles 
correspond to the correlation function between haloes and neighbors in the direction perpendicular to the angular momentum of the 
centre haloes. The diamonds show the results when using tracers along the direction of the angular momentum, and the squares the 
results when all neighbors are taken into account. The lower left panel shows the ratio between the correlation functions shown in the 
top-left panel and the results using all neighbors. Right panel: same as left panel, for a high mass DM halo subsamplc from the large 
simulation. 



triangles and squares, respectively). The lower panel shows 
the ratios between these three estimates and the correlation 
function considering all neighbors; their relative amplitudes 
help identify alignments in either direction. As can be seen, 
the relative amplitudes of the correlation function when us- 
ing neighbors in the directions parallel and perpendicular 
to the angular momentum of the DM halo are different; for 
small separations R < lh~^Mpc, there is an excess of cor- 
relation in the direction perpendicular to the angular mo- 
mentum of up to a 20% effect (between correlations parallel 
and perpendicular to the angular momentum, lower panel), 
whereas for larger separations the alignment diminishes and 
shows a tendency to reverse (although, as can be seen in 
the figure, errors are larger at these separations); we come 
back to this later in this section where we use a higher signal 
estimator for alignments. 

When using higher mass DM haloes, the effect is also 
present with an even higher difference and statistical signif- 
icance between the structure in the direction of the angular 
momentum and perpendicular to it (Right panels of Figure 
[l]). In this case, neighbors in the direction perpendicular to 
the angular momentum show more correlation, and there is 
a clear tendency for this alignment to decrease with separa- 
tion, which goes from almost a 50% effect on relative ampli- 
tudes (lower panel) at small separations to about < 5% at 
large separations. Note that the alignment does not change 
at large separations as seemed to be the case for low mass 
haloes, even at _R = 20h~^Mpc compared to R = 5h~^Mpc 
analised in the low mass halo sample. 

Notice that the halo-particle cross-correlation functions 
in Figure [T] show the characteristic 1- and 2-halo terms re- 



ported in several previous measurements of both observa- 
tional (e.g. Zehavi et al., 2004, Cooray, 2005) and numerical 
simulations (Zheng et al., 2004); these two regimes are sep- 
arated in the figure by the vertical dashed lines. This sepa- 
ration is obtained by measuring the scale at which the cor- 
relation function is seen to shift between these two regimes 
(there is a local minimum in the first derivative of the cor- 
relation function). It can be seen that the alignment signal 
shows a clear transition from the 1- to the 2-halo regime; in 
particular, low mass haloes show indications of a different 
angular momentum— large-scale alignment behaviour in the 
2-halo regime, suggesting an important relationship between 
alignment and halo term. This can be interpreted as differ- 
ent physical mechanisms producing the alignment in the 1- 
and 2-halo terms. 

We make quantitative estimates of the alignments for 
the 1- and 2-halo terms separately from the cross-correlation 
functions shown in Figure [T] as well as from the other sam- 
ples of DM haloes with intermediate masses (for a total of 
13 samples of varying halo masses). We only consider pairs 
that are separated by a distance at least as large as the tran- 
sition between regimes for the 2-halo term, and at shorter 
separations for the 1-halo calculations. 

In order to avoid effects of covariance between the cor- 
relation function bins to some extent, we make global 1- 
and 2-halo term estimates of the ratios between the correla- 
tion functions in the directions perpendicular and parallel to 
the angular momentum of haloes, using all the halo-particle 
pairs separated by distances within the 1- and 2-halo ranges 
(i.e. we do not use the narrow bins in log]^g(r) used to calcu- 
late the cross- correlation functions). We acknowledge that 
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traces of covariance will still be present between the 1- and 
2-halo term ratios. We present these results in Figure [2j 
where triangles show the 1-halo term ratio between correla- 
tion functions for neighbors in the directions perpendicular 
and parallel to the halo angular momentum; Diamonds cor- 
respond to the 2-halo term. The horizontal long dashed line 
C||/^-L = 1 corresponds to the isotropic case. The data- 
points at masses lower than 10^^ Mq correspond to the 
small numerical simulation, and show larger errors due to 
both, the smaller simulated volume, and the lower ampli- 
tude of clustering characterising low mass haloes. 

As can be seen, the 1-halo term always shows more 
structure in the direction perpendicular to the angular mo- 
mentum, and more so for larger masses, whereas the 2-halo 
term shows a similar behaviour with the possible exception 
of low mass haloes, whose angular momenta tend to point 
in the direction of the large-scale structure (although with 
a low statistical significance). The alignments are detected 
at about a 5 — a level for the higher mass halo subsamples 
drawn from the large simulation. As in the 2-halo regime, 
low mass haloes show considerably higher errors. 

In the case of low mass halos {M ~ 10^'^/i~^Mq) in 
the 2-halo term regime (diamonds in Figure [2]), the align- 
ment signal is the opposite to that of higher mass haloes (a 
10 ± 5% effect). This anti-alignment signal fits in well with 
the overall tendency for this quotient S.\\/(,± to increase with 
mass, which can be fit by a log-linear relation which crosses 
the unit value at ~ 10^^ Mq. This mass limit suggest 
a change in the alignment signal in qualitative agreement 
(regardi ng the alignment di rection) with results on simula- 
tions bv iHahn et al.l (j2007al ). with masses above and below 
the typical mass scale for collapse M* = 5.5 x 10^'^'^ Mq 
at z = 0. The best-fit log-linear relation for the 2-halo term 
results is shown in Figure [2] (dash dotted line), and corre- 
sponds to, 

^||/C± = (0.95±0.01)-h(0.082±0.008)*[/ogio(M)-11.5] (2) 

For the 1-halo term, the dependence of alignment on mass 
is well adjusted by (dashed line on Figure [2} , 

Cll/Cx = (0.95 ± 0.01) (0.082 ± 0.008) *[;o5io(Af)- 11.5] (3) 

The 2-halo term effect can be understood within the 
framework of TTT, specially for the high mass haloes, which 
readily show that their angular momentum is perpendicular 
to the present-day mass distribution. It seems that for high 
enough mass haloes, the local large-scale structure still de- 
scribes the original tidal field. In the case of low mass haloes, 
the alignment seems to reverse. This result can be under- 
stood by an intuitive non-linear effect beyond the abilities 
of TTT, for example if these haloes have formed recently in 
filaments, from the merger of material arriving from voids 
surrounding these elongated structures. 

The 1-halo term results indicate that regardless of halo 
mass, the structure within dark-matter haloes is preferen- 
tially aligned with the plane perpendicular to the angular 
momentum, with a higher statistical significance and signal 
than the 2-halo term. 
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Figure 4. Top panel: Correlation function around "edge-on" 
red spiral galaxies in the SDSS-DR6 (sample 03), in the direc- 
tions parallel and perpendicular to the line of sight. The colour 
indicates the amplitude of the correlation function (see the fig- 
ure key). Bottom panel: real-space correlation function for the 
same sample of galaxies, obtained as explained in the text (sym- 
bols with error bars) , and the best- fit ACDM real-space correlation 
function (black lines). The dashed lines show the range of corre- 
lation function amplitudes allowed by the errors in the data. The 
gray thick line shows the mass real-space correlation function in 
the ACDM model. 

2.2 Simulating observational biases in the 
numerical simulation 

It can be difficult to estimate the angular momentum of real 
galaxies. In the case of the SDSS-DR6, which we use in this 
work, the only information that can allow us to obtain an 
estimate of the direction of the angular momentum is the 
photometry of galaxy images, specially for spiral galaxies 
but also for ellipticals. We assume that spiral discs or fiat- 
tened spheroids are perpendicular to the angular momentum 
of galaxies and therefore, edge-on galaxies (either disks or 
spheroidals) will be more likely to have their angular mo- 
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Figure 3. Left panels: projected correlation functions around "edge-on" low mass haloes in the small numerical simulation (top-left), 
in the directions perpendicular and parallel to the angular momentum of centre haloes (triangles and diamonds symbols, respectively). 
Ratios between these projected correlation functions are shown in the bottom-loft panel. Right panels: same as left panels for high mass 
haloes in the large numerical simulation. 



mentum on the plane of the sky. Then, one can calculate 
the cross-correlation function for neighbors in the projected 
direction parallel and perpendicular to the direction of the 
angular momentum of the centre galaxies, in redshift slices 
around them, assuming all these neighbors to be at the same 
distance than the centre galaxy. This slice needs to be de- 
fined so that it avoids projecting large amounts of foreground 
and background structures while still allowing enough num- 
bers of neighbors to obtain good statistics. 

In order to assess whether the alignments found in the 
previous subsection are identifiable in observational cata- 
logues analised this way, we reproduce the procedure men- 
tioned in the previous paragraph in our analysis of the nu- 
merical simulation which consists on projecting the struc- 
ture onto the plane of the sky, assuming all objects within 
An = 750km/s from the centre haloes to be exactly at the 
same distance from the observer. 

Therefore we only use information from two cartesian 
coordinates in the simulation, x and y, and transform the 
third, z, into a velocity component by multiplying z by the 
Hubble Constant and adding the z component of the pecu- 
liar velocity. In order to complete the emulation of redshift 
space effects suffered by observational surveys, we apply this 
transformation to the simulation particles as well as to the 
halo catalogue. We then select "edge-on" haloes from the 
simulation, by requiring that their angular momentum is 
at an angle less than 60° from the sky, represented by the 
X — y plane 0- We can then study two projected correlation 
functions centred in these objects using DM particles, (i) 



^ Such angle threshold is equivalent to a ratio of 0.5 of the ellip- 
soidal image semi-axes for a perfect thin disc. 



a correlation function calculated using pairs that subtend 
an angle with the projected angular momentum of less than 
30°, which we call the parallel projected correlation function 
(tj||), and (ii) a perpendicular projected correlation func- 
tion {u)_l) obtained by counting pairs that subtend an angle 
greater than 60° from the projected angular momentum. 

Figure |3] shows the results from this procedure. The left 
panels show the results for low mass haloes in the small sim- 
ulation (projected correlation functions in the top- left, and 
ratios in the bottom- left), whereas the right panels show 
the results for high mass haloes in the large simulation. As 
can be seen, these results are qualitatively consistent with 
the alignment signal found in the previous subsection using 
the full three-dimensional correlation function, which indi- 
cates that projected correlation functions are a suitable tool 
to detect this effect. The overall effect has a lower ampli- 
tude; high mass haloes show a ratio between structure in 
the directions perpendicular and parallel to their angular 
momenta of about a 30% excess (for the 3D case, this was 
above 50%) and low mass haloes show a 10% effect (com- 
pared to 20 — 30%), which makes this detection a slightly 
more difficult task. 



3 OBSERVATIONS 

Galaxies with spiral and flattened spheroidal shapes pro- 
jected onto the sky give a handle on their orientation with 
respect to the line of sight and, therefore, it is possible to 
estimate in a statistical sense the direction of their angular 
momentum. Spiral galaxies with round projected shapes are 
likely discs seen face-on and therefore their angular momenta 
point in the direction perpendicular to the plane of the sky; 
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Table 1. Observational samples from the SDSS-DR6 spec- 
troscopic main galaxy sample. The first and second column 
indicate the sample name and number of galaxies in them, 
the third column indicates the number density, and the fourth 
column the host-halo mass. 



Sample 


Members 


n/lO-^h-^Mpc^ 


Logio(Mhost/h~^Af0) 


Ol 
02 
03 
04 
05 


132000 
59300 
72700 
10900 
1400 


10.4 ±0.8 
6.50 ±0.12 
5.50 ±0.12 
7.15 ±0.17 
1.10 ±0.20 


13.06+o« 

ii.^U_Q gg 

12.19l?f7 



elongated projected shapes indicate discs seen edge-on, and 
therefore their angular momenta point in the direction per- 
pendicular to the long axis of the ellipsoid. On the other 
hand, elliptical galaxies with round projected shapes can ei- 
ther be nearly intrinsically spherical or flattened spheroids 
seen face-on; however, elongated ellipticals are more likely 
flattened spheroids seen edge-on since recent results on in- 
trinsic elliptical shapes indicate that these are more com- 
monly oblate rather than prolate spheroids (e.g. Padilla & 
Strauss, 2008). Therefore, elongated ellipticals are likely to 
be also characterised by angular momenta pointing in the di- 
rection of the plane of the sky roughly perpendicular to the 
major axis of the ellipsoid, as is the case for edge-on spiral 
galaxies. These last two types of galaxies will be the centre 
of the following analysis. We notice that the amplitude of 
the angular momentum of flattened spheroidals is likely to 
have a low amplitude relative to that of spiral galaxies. 

The sample we selected to carry out our analysis is 
the spectroscopic SDSS-DR6, which contains a total of 
~ 580, 000 galaxies with spectroscopic redshifts and pho- 
tometry in five bands, u,g,r,i, z, as well as parameters de- 
termining the projected galaxy shapes consisting of the two 
semi-axes and the position angle of the ellipse that once 
convolved with the PSF (seeing) provides the best match to 
the photometric image of each individual galaxy. In the re- 
mainder of this paper, we only consider SDSS-DR6 galaxies 
within the redshift range 0.02 < z < 0.09. 

We study the alignment signal in the observational data 
for different subsamples defined by galaxy colour and lumi- 
nosity in order to detect trends that can be associated to the 
results from the numerical simulation. The subsamples we 
study are sample 01, consisting of all the edge-on galaxies 
with axial ratio b/a < 0.7 in the SDSS-DR6; 02, red galax- 
ies with g — r > 0.7; 03, blue galaxies with g ~ r < 0.7; 04, 
faint galaxies with Air > —19.5; and 05 composed by bright 
galaxies with A'lr < —19.5. The selection of subsamples is 
done so as to ensure a number of target galaxies to allow 
good statistics (> 1000); the number of galaxies in each 
sample is shown in Table [1] where we also show the galaxy 
number density calculated out to the redshift of complete- 
ness of each sample. 

In order to characterise in a quantitative way the sub- 
samples selected from the SDSS-DR6, we calculate the me- 
dian mass of ACDM haloes predicted to show the same 2- 
halo regime clustering amplitude as the galaxies in each sub- 
sample. We follow the procedure used in several previous 
works (cf. Croft et al. 1999, Padilla et al., 2001) consist- 
ing on calculating the auto-correlation function of galaxies 



in two coordinates, one parallel {a) and one perpendicular 
(tt) to the line of sight, f((T, tt), and integrating over the tt 
direction to avoid the effect of redshift space distortions. 
The result from this integration is the projected correlation 
function (see for instance. Croft et al., 2001), H(cr), 

E{a) = 2 [ "^^^ ^{a,Tv)d7v, (4) 

which in turn can be inverted to obtain the real-space cor- 
relation function via 

where the sum is performed over the bins in a where the pro- 
jected correlation function has been calculated. The masses 
of host haloes can then be obtained by finding the bias be- 
tween the resulting real-space correlation functions and the 
matter correlation functions in a ACDM cosmology (both in 
the 2-halo regime), which we obtain by Fourier transform- 
ing the non-linear matter power spectrum (from Peacock & 
Smith, 2000) for the same cosmological parameters used in 
the numerical simulations. With an estimate of the bias, the 
mass of the host haloes can be found by using the Sheth, 
Mo & Tormen (2001) formalism. Note that this procedure 
is equivalent to measure the real-space halo-halo correlation 
function for haloes of different masses (Sheth, Mo & Tor- 
men, 2001, Padilla & Baugh, 2002), with the advantage of 
avoiding large uncertainties from small samples of haloes re- 
stricted to narrow ranges of halo mass in the simulation. We 
show one example of this procedure in Figure |4l where we 
find the mass of haloes with an equivalent clustering ampli- 
tude to the sample of red galaxies in the SDSS-DR6 (02). 
In this figure, the top panel shows the auto-correlation func- 
tion, ^((T, tt), of 02 galaxies, and the lower panel the result- 
ing real-space correlation function (symbols) and the best 
fit ACDM curve. From this analysis we find that red galax- 
ies in the SDSS-DR6 (with b/a < 0.7) are consistent with 
haloes of a median mass of log-^^Q{M/h~^ Mq) = 13.06lg 27- 
This provides us with a quantitative way to compare the 
SDSS-DR6 subsamples to the halo samples in the numerical 
simulations. Also, in order to infer this scale for our obser- 
vational samples, we will use the dependence of the 1- to 
2-halo term transition on DM halo mass measured in our 
the numerical simulations. 

The other SDSS-DR6 galaxy samples are fitted by the 
correlation function of DM haloes of masses indicated in 
Table [T] Notice that our study covers over one and a half 
decades in the host-halo mass of SDSS-DR6 galaxies. 

As we demonstrated in the previous section both, the 
large scale structure and alignments surrounding galaxies, 
can be studied using the projected correlation function, 
io{a), measured in terms of the projected separation, a. We 
now apply this method to our samples of edge-on galaxies, 
and compare the outcome from using tracers in the direction 
parallel and perpendicular to the inferred angular momen- 
tum, to detect differences between the structure along these 
two directions. 

The projected correlation functions are calculated us- 
ing galaxies in a given subsample, cross-correlated to the full 
spectroscopic SDSS-DR6 catalogue. As in the 2-dimensional 
analysis carried out in the numerical simulations, we cal- 
culate the cross- correlation for edge-on galajcies selected so 
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Figure 5. Correlation functions around "edge-on" galaxies in 
the SDSS-DR6, in the direction perpendicular and parallel to 
the angular momentum of centre haloes (triangles and diamonds, 
respectively). Top panel corresponds to sample Ol (full SDSS- 
DR6), middle panel to 02 (red galaxies), and bottom panel to 
sample 03 (blue galaxies). The vertical dashed lines show the 
transition from the 1- to the 2-halo regimes corresponding to the 
host-halo mass obtained from the real-space correlation function 
of the different samples. 



that h/a < 0.7. Note that in the numerical simulation we 
used an angle between the halo angular momentum and the 
line of sight of 60°, which in the case of a perfect thin disk 
is equivalent to b/a = 0.5, but is larger for a thick disk as is 
the case of the intrinsic shapes of spiral galaxies or flattened 
ellipsoids (Padilla & Strauss, 2008). 



Figure [5] shows the measured correlation functions for 
Edge-on galaxies for the three subsamples that show signif- 
icantly different signatures of alignments with the angular 
momentum. These are the full sample of edge-on spirals, 
01 (top panel), the red galaxy sample, 02 (middle panel), 
and the blue galaxies, sample 03 (bottom panel). As in the 
results for the numerical simulation, triangles show the cor- 
relation with galaxies in the direction of the angular momen- 
tum, diamonds in the perpendicular direction, and squares, 
when using all neighbors. A first glance at the shapes of 
these correlation functions indicates that there is not as clear 
a transition between the 1- and 2-halo terms as in the re- 
sults from the numerical simulations (see Fig. [Sjl; in these 
samples the main effect is a transition to a roughly constant 
slope for the log(5) vs. log a relation. This may be an indica- 
tion that the internal structure of haloes is not traced in the 
same way by galaxies in the SDSS-DR6 and DM particles 
in the ACDM model. Regarding the alignments, it can be 
seen that the three subsamples show the large-scale struc- 
ture preferentially aligned with the direction perpendicular 
to the angular momentum, specially for the 2-halo regime, as 
was found in the cosmological simulation. The difference be- 
tween the two correlation functions is roughly a 2 — cr detec- 
tion over the range of scales defined by 2 < r/h~^Mpc< 20 
and corresponds to about a 2 — 4% effect. 

The analysis of the faint and bright samples (04 and 
05, respectively) provide alignment results in quantitative 
agreement with the full galaxy sample, 01, with results 
within a 1 — a difference. Therefore, we will not show re- 
sults from these two samples in the remainder of this work. 

A more quantitative comparison between observational 
and simulation results is presented in the following section. 



4 DISCUSSION 

In order to compare quantitatively the alignments found in 
the observations and in the numerical simulation, we use the 
transition scale from the 1- to the 2-halo regimes inferred in 
the previous section to determine the degree of alignment 
found in these two regimes separately. As in section [2711 we 
estimate ratios between the correlation functions in the di- 
rections perpendicular and parallel to the estimated galaxy 
angular momentum direction, using all pairs within the 1- 
and 2-halo regimes. This can then be directly compared to 
the results from the numerical simulations for projected cor- 
relations. 

Figure shows the ratios between projected correla- 
tion functions in the directions parallel and perpendicular 
to the angular momentum as a function of halo mass. The 
results from the numerical simulations are shown as dia- 
monds, and the results from the observational samples at 
the corresponding host-halo masses are shown as filled sym- 
bols. We show results for samples 01 (full SDSS-DR6, filled 
squares), 02 (red galaxies, filled upward pointing triangles), 
and 03 (blue galaxies, filled downward pointing triangles). 
The top panel shows the results from the analysis of the 

1- halo term, the lower panel from the 2-halo term. For the 

2- halo term regime the resulting alignment signals for ob- 
servational samples are 1.8 ± 0.5% for the full sample Ol, 
2.7±0.7% for the red galaxy sample 02, and 0.8±0.4% for the 
blue galaxy sample 03. The 1-halo term regime shows align- 
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Figure 6. Ratios between projected correlation functions in the 
directions parallel and perpendicular to the projected angular 
momentum as a function of halo mass in the numerical simu- 
lation (open diamonds). Galaxy samples Ol (all galaxies), 02 
(red galaxies), and 03 (blue galaxies), are shown as filled sym- 
bols: squares, upward pointing triangles, and downward pointing 
triangles, respectively. Top panel: ratios for the 1-halo term. Bot- 
tom panel: ratios for the 2-halo term. 



ments of 0.1 ± 0.5%, 0.1 ± 0.5%, and 1.0 ± 0.4% for samples 
Ol, 02, and 03, respectively. As can be seen, the alignment 
found in the 2-halo term is consistent between observational 
samples and the numerical simulation (differences are at a 
maximum a 1 — ct effect). 

This agreement indicates that TTT in combination with 
hierarchical clustering are a viable theory since not only the 
angular momentum of galaxies points in the direction per- 
pendicular to the surrounding large-scale structure, but also 
the amplitude of the effect is similar to what is found in the 
numerical simulations. It should be noticed that this would 
also indicate that the direction of the angular momentum 



in galaxies in these samples is tracing the angular momen- 
tum of the host halo. The results for the 1-halo term, on the 
other hand, show that galaxies follow a different behaviour 
than DM haloes as was suspected from the lack of a clear 
transition to the 1-halo term in the shape of the SDSS pro- 
jected correlation functions; at the range of masses corre- 
sponding to the observational samples, the simulated haloes 
show very significant alignments, whereas the only sample of 
galaxies that a shows non-zero alignment is that composed 
by blue galaxies (although still lower than the simulation 
data, at a regime of low expected alignments). The combi- 
nation of the 1- and 2-halo term results could indicate that 
the general galaxy population in the SDSS-DR6 which our 
samples are cross-correlated to, can be used as good tracers 
of the large-scale structure (i.e. have not suffered important 
changes to deviate from the angular momentum of their host 
haloes) but would fail to serve as indicators of the internal 
distribution of matter in their host DM haloes; different as- 
trophysical processes could produce this effect, specially for 
such low M ~ W^^ Mq halo masse^. This result is par- 
ticularly interesting since blue galaxies do show traces of an 
alignment in the same direction as the numerical simulation, 
which has a number of possible explanations. For instance, 
the angular momentum of red galaxies could be subject to 
a reshuffle that adds enough noise to loose its alignment 
with the internal halo structure while leaving the large-scale 
signal almost intact, an effect could also involve a positional 
reshuffle within the halo structure; another possibility is that 
red galaxies are able to rearrange their companion galaxies 
inside DM haloes. However, a more quantitative analysis 
is required in order to find the true underlying cause for 
this behaviour. We acknowledge the possibility that obser- 
vational biases not included in our analysis of the numerical 
simulation are still affecting these measurements. 



5 CONCLUSIONS 

We study the alignments between the angular momentum of 
individual objects and the large-scale structure in cosmolog- 
ical numerical simulations and in the SDSS-DR6. The angu- 
lar momentum of DM haloes in ACDM simulations is found 
to be preferentially oriented in the direction perpendicular 
to the distribution of matter in the range of separations cor- 
responding to the 1- and 2-halo terms. These results are in 
agreement with the Tidal Torque Theory. 

We find that more massive haloes show a higher de- 
gree of alignment, and that the 1-halo term shows a much 
higher alignment than the 2-halo term. We find a maximum 
1-halo alignment of ~ 40% for haloes of ~ lO^^'/i'^M©; at 
low masses, the 1-halo signal seems to tend to a constant 
value of 20%. The 2-halo term alignment shows a maximum 
signal of ~ 15% at high masses, and possible inversion at low 
masses where the angular momentum might even preferen- 
tially point in the direction of the distribution of mass on 
large scales. The latter would be consistent with low mass 

There is observational evidence that groups of galaxies with 
masses above this limit are characterised by intr insic shapes con - 
sistent to those found in numerical simulations jPaz et al]|2006h : 
however it is difficult to directly measure the shapes of lower mass 
systems due to discreteness effects. 
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groups formed not too long ago from mass distributed in the 
surrounding areas of the already-formed filaments or walls. 

The alignment of the angular momentum with the 
large-scale structure found for l ow mass systems is in agree- 



ment with previous works by (jTruiillo. Carretero fc Patiri 
l2006h o n redshift su rveys, a nd by iBrunino et al 
I2OO3); iPatiri et all (|2006[): ICuesta et alj j2008h 



lAragon-Calvo et al.l (12007): Hahn et all (|2007al l 



simulations. The main difference with these previous stud- 
ies relies in that we quantify the strength of the alignment 
in terms of the two-point correlation function, as a function 
of halo mass. In particular, we find that a log-linear relation 
fits the dependence of alignment on mass, for both, the 1- 
and 2-halo term regimes. For the latter we observe a turning 
point between an angular momentum pointing in the direc- 
tion perpendicular to the structure, to no-alignment, and 
then to a regime where the angular momentum is contained 
by the plane defined by the surrounding structure. This 
occurs at a mass ~ M«, in qualitat i ve agr e ement with 
previo us studies (jBailin fc Steinmet j l2005l : iHahn et al.l 
l2007al lb[). 

Combining our results for the 1- and 2-halo terms for 
haloes in the numerical simulation, namely that the angular 
momentum of haloes is perpendicular to the bulk of the mat- 
ter within haloes (1-halo term), and also to the surrounding 
large-scale structure (for masses > M"), our findings are 
in agreement with previous works who find that the major 
axis of haloes is usually aligned with the large-scale struc- 
ture (Gottlober & Turchaninov 2006, Rett et al. 2007, Fal- 
tenbacher et al. 2002; Kasun & Evrard 2005; Bailin & Stein- 
metz 2005; Colberg, Krughoff & Connolly 2005; Altay, Cold- 
berg & Croft 2006; Basilakos et al. 2006; Ragone-Figueroa & 
Plionis 2007, and references therein). Regarding the 1-halo 
term signal, our resu lts are also i n q ualitative agreement 
with recent results by Kne be et al.l |2008) on the alignments 
of substructure within dark matter halos. 

We use numerical simulations to reproduce the proce- 
dure that can be applied to observational samples using pro- 
jected correlation functions. We find that the alignment sig- 
nal diminishes considerably but that it can still be measured. 

Therefore, we study edge-on galaxies in the SDSS-DR6, 
and assign the direction perpendicular to the major semi- 
axis as the direction of the angular momentum. Since our 
samples are restricted to objects with fiattened apparent 
shapes {h/a < 0.7), these include both edge-on spiral galax- 
ies and flattened spheroids. We are able to detect alignments 
in all our galaxy samples. These include the full SDSS-DR6 
edge-on galaxies, a sample composed by blue {g — r < 0.7) 
galaxies, the sample of red {g — r> 0.7) galaxies, and a sam- 
ple of faint and another composed by bright galaxies. The 
latter two samples show a quantitative agreement with the 
results from the full sample of edge-on galaxies in the SDSS. 

In all cases we find a significant excess of structure in 
the direction perpendicular to the angular momentum for 
the 2-halo term, in good agreement with the numerical sim- 
ulation results. The 1-halo term results show either a null 
or lower alignment signal than expected from the analysis 
of projected correlations in the numerical simulation. This 
would suggest the effects of astrophysical processes acting so 
that galaxies do not follow the DM structure inside low mass 
haloes of W^^h~^ Mq or lower masses, a result that com- 
plements previous studies of high mass groups in the SDSS 



with M > 10^"^ Mq, where galaxies are found to follow 
the internal DM structure (e.g. Paz et al., 2006). This effect 
is apparently more important for red, g—r > 0.7 galaxies, in- 
dicating that there may be interesting astrophysical effects 
at work at low mass dark-matter haloes which could pro- 
duce changes in the orientation of galaxy angular momenta, 
or the spatial re-distribution of galaxies within them. 
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